Density Fluctuation Spectrum of Solar Wind Turbulence between Ion and Electron Scales 
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We present a measurement of the spectral index of density fluctuations between ion and electron scales in 
solar wind turbulence using the EFI instrument on the ARTEMIS spacecraft. The mean spectral index at 1 AU 
was found to be -2.75 ± 0.06, steeper than predictions for pure whistler or kinetic Alfven wave turbulence, but 
consistent with previous magnetic field measurements. The steep spectra are also consistent with expectations of 
increased intermittency or damping of some of the turbulent energy over this range of scales. Neither the spectral 
index nor the flattening of the density spectra before ion scales were found to depend on the proximity to the 
pressure anisotropy instability thresholds, suggesting that they are features inherent to the turbulent cascade. 

PACS numbers: 94.05.Lk, 52.35.Ra, 96.60.Vg, 96.50.Bh 
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Introduction. — The solar wind has been observed to be a 
turbulent plasma for many decades (see |l]-|5| for recent re- 
views). Understanding its properties is important for deter- 
mining the universal features of turbulence and how the solar 
wind and collisionless plasmas in general are heated 0, Q] . 
This Letter examines the properties of small scale density fluc- 
tuations in the solar wind at plasma kinetic scales, where dis- 
sipation of the turbulent energy is thought to take place. 

In the inertial range, i.e., scales larger than the ion ki- 
netic scales, the one-dimesional magnetic field power spec- 
trum P(k) ~ k a , where k is the wavenumber, is observed 
to have a spectral index close to a = -5/3 (e.g., ]8[]), con- 
sistent with a turbulent cascade. It has been known for many 
years 10-12] that this spectrum steepens around ion kinetic 
scales, although it is still not yet obvious at which ion scale 
the steepening occurs lfl3ll . More recently, a further change in 
the spectrum has been reported around electron scales, with 
either a steeper power law lUill or exponential falloff 1U5I sug- 
gested. Early measurements of the spectral index between ion 
and electron scales showed a wide range between -4 and -1 
although more recently values between -2.9 
and -2.3 have been obtained iTfl [HI [H [HI . 

The steepening at ion scales was originally attributed to ion 
cyclotron damping [12, T^, 2^1 but it was later suggested that 
the dispersive nature of fluctuations at these scales could also 
be the cause 121U23I1 . It was proposed that the power law be- 
tween ion and electron scales could be explained by a turbu- 
lent cascade mediated by the dispersive fluctuations 124142611 
similarly to the Alfvenic cascade at larger scales []27|-|29J . 

Theoretical predictions of t he sp ectral index of the disper- 
sive cascade have been mad e 125 based on Kol- 
mogorov scaling arguments foil . If the turbulence is strong 
(non-linear eddy timescales « linear wave timescales), the 
magnetic field spectral index is predicted to be -7/3. Addi- 
tional effects can be included to account for the steeper ob- 
served spectra, for example, shear generated cyclotron reso- 
nant waves [34], an ion entropy cascade 126I1 . wave-particle 
scattering 13511. electron Landau damping [36], nonlocal inter- 
actions 13711 or increased intermittency 03811 . 

Since more than one type of plasma wave can exist at these 
scales, the nature of the dispersive cascade is debated. It is 
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kinetic Alfven waves (KAWs) UJ, [1J, |26L |3(| [H ErH. 
Since both wave modes produce turbulence with the same 
spectral index (-7/3), other tests have been used to distinguish 
between them 1 42 , 46|- 49(] . Other possible contributions to the 
spectrum at these scales include current sheets 0111 and 
kinetic instabilities 15211 and their effect remains to be fully 
investigated. 

The spectrum of density fluctuations has been well mea- 
sured in the inertial range (e.g., 1531, 15411 ) but since current par- 
ticle counting instruments take a several seconds to generate 
a density moment, it is not currently possible to measure the 
density spectrum below ion scales with this technique. Higher 
frequency measurements from the ISEE propagation experi- 
ment 15511 show the density spectrum flattening before the ion 
scales then steepening at smaller scales, although the steep- 
ening was attributed to the measurement technique. Similar 
spectra were seen with Cluster [1561 1 using the spacecraft poten- 
tial measurement as a proxy for density (as described below), 
although the data resolution was not sufficient to measure far 
beyond the ion scales. Radio scintillation measurements also 
suggest a steepening of the density spectrum at ion scales in 
the inner solar wind [57] and interstellar medium [58 1. 

In this Letter, we present new measurements of the solar 
wind density spectrum at 1 AU that have a low noise level and 
sufficient resolution to allow the spectral index between ion 
and electron scales to be determined. 

Measurement technique. — In sunlight, spacecraft emit pho- 
toelectrons and typically become positively charged. This 
attracts a return current of electrons from the surround- 
ing plasma, reducing the spacecraft potential relative to the 
plasma, V sc , until an equilibrium is reached in which the cur- 
rents to and from the spacecraft are balanced. For higher 
electron density, n e , the return current is larger, resulting in 
smaller U sc . Thus, V sc can be used as a proxy for n e 15911 . 
allowing density fluctuations to be measured at a higher fre- 
quency than with particle counting instruments. 

V sc is a good proxy for density at frequencies lower than 
the inverse time it takes the spacecraft to charge, which is 
determined by dV sc /dt = I t /C, where C is the spacecraft 
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capacitance and I t is the total current to the spacecraft. The 
important contributions to I t are the photoelectron current and 
plasma return current, 7 p i, giving/ t w 7 p ] — 7 pe e _1/sc /' r P c , where 
Ipe is the photoelectron current at = and T pe is the pho- 
toelectron e-folding energy. Applying small perturbations to 
the equilibrium I p \ and V^, it can be shown that the spacecraft 
relaxes exponentially to the new equilibrium in response to 
density changes with time constant « CT pe /I p \. This corre- 
sponds to a frequency w 6 kHz in the solar wind and in this 
Letter, we consider fluctuations at much lower frequencies, 
where V^ c can be well calibrated to n e . 

Several intervals of high frequency data from the 
ARTEMIS -P2 spacecraft [60j] were used, which were in the 
free solar wind [161fl and for which a reliable conversion from 
V sc to n e could be made. V^ c is measured by the EFI instru- 
ment Jill], which consists, in part, of four conducting spheres 
coupled to the plasma at the end of orthogonal booms in the 
spacecraft spin plane. In this Letter, data from one pair of op- 
posite probes (probes 1 and 2) was used, since the others were 
found to contain large spin period spikes in their time series, 
likely due to shadows from one of the axial booms momentar- 
ily altering the probes' photoemission. The measured poten- 
tial of the two probes relative to the spacecraft was averaged 
to reduce offsets due to solar wind electric fields. 

The probes themselves also charge positive from photoe- 
mission and are supplied a bias current to reduce their po- 
tential, but are left to remain about 1 V higher than the sur- 
rounding plasma. This places them at a point on their current- 
voltage curve where their potential is far less sensitive to den- 
sity fluctuations than the spacecraft potential is. In addition 
to this 1 V offset, a further scale factor correction of 1 . 15 was 
applied to convert the average probe potential to a measure- 
ment of V" sc . This accounts for the fact that the probes are not 
infinitely far from the spacecraft but measure plasma which is 
slightly perturbed by the spacecraft environment (see Section 
2.1 of 116311 for details of these corrections). 

To obtain a calibration curve to convert from Vx to n e , 
spin resolution K c data was compared to n e data from the 
ESA instrument l63ll . An example of this comparison for 
the interval 00:04 - 02:30 on 1 1th October 2010 is shown in 
Fig.Q] The electron density from ESA was estimated from the 
measured ion density, assuming that 4% of the ions were al- 
phas and the rest protons (see Section 3.2.1 of JmI). Since 
n e is expected to be proportional to the exponential of V sc 
II591 16511 . a least squares fit of the data in Fig.[T]to the equation 
n e = Pi cxp ({Vsc — P2) /P3) was performed, where Pi, P2 
and P3 are fit parameters. The P2 parameter was included to 
allow for variations in the 1 V probe potential offset. The fit 
is shown as a red line in Fig.Q] 

Density spectrum. — The calibration curve from Fig. [T]was 
applied to the 128 samples/s data obtained during the 
"particle burst" mode interval 00:21 - 01:14 on 11th Octo- 
ber 2010 to obtain a density time series n e (t). The power 
spectrum of density fluctuations as a function of spacecraft- 
frame frequency is given by P(/ sc ) = /_ R(T)e~ 2m * acT dr , 
which is the Fourier transform of the autocorrelation func- 
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FIG. 1. Electron density, n e , as a function of the spacecraft potential 
relative to the solar wind, V sc . The best fit exponential calibration 
curve is shown in red. 



tion R(t) = (n s (t)n e (t + r)), where the angular brackets 
denote an ensemble average. The power spectrum was es- 
timated using the multitaper technique with time-bandwidth 
product NW = 4 ll66ll and is shown in Fig. [2^ in blue. In the 
same figure, the spectrum of a 8192 samples/s "wave burst" 
mode interval 00:36:01 - 00:36:05 on 11th October 2010 is 
shown in green. 

Several features can be seen in the spectrum. Large spikes 
at harmonics of the spacecraft spin frequency (0.30 Hz) are 
present throughout the spectrum. These are caused mainly by 
the varying illumination of the grounded sections of the EFI 
booms as the spacecraft spins, altering the spacecraft photoe- 
mission and, therefore, the spacecraft potential ]62|]. All in- 
tervals in this Letter were reduced to an integer number of 
spin periods to reduce spectral leakage from the spin harmon- 
ics. Since they are relatively localized in frequency, these 
harmonics were eliminated by removing up to 0.03 Hz either 
side of the spin harmonics, plus two additional sections: 1.07 
Hz - 2.28 Hz and 1.35 Hz - 2.52 Hz. Similarly, spikes in 
the spectrum at harmonics of onboard clocks operating at 8 
Hz and 32 Hz 16711 were removed. The resulting spectrum is 
shown in Fig. |2j?. The low frequency portion of the high reso- 
lution (green) spectrum has also been removed since leakage 
from the spin harmonics here is large due to the short interval 
length. The spectrum was also smoothed by averaging in 45 
logarithmically spaced bins from 2 x 10~ 3 Hz to 1 x 10 3 Hz 
(Fig. |S). 

Under Taylor's hypothesis ]68|], the measured frequency 
spectrum can be interpreted as a wavenumber spectrum since 
the spacecraft-frame frequency is / sc = fcv sw /(27r). This re- 
quires the fluctuation speeds to be less than the solar wind 
speed, which is well satisfied for Alfvenic turbulence in the 
inertial range but may or may not be valid below ion scales. 
There is mounting evidence from phase speed 11421 14811 and 
polarization i47ll measurements that the fluctuations between 
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FIG. 2. Power spectra of electron density fluctuations: (a) from cal- 
ibrated data (artificial spikes present), (b) with artificial spikes re- 
moved, (c) smoothed. Ion and electron gyroscales, p, and inertial 
lengths, d, are marked with vertical dashed lines. The empirical noise 
floor is marked with a horizontal dotted line. 



ion and electron scales are KAW-like. Since KAWs are low 
frequency (compared to the ion cyclotron frequency), this sug- 
gests a wavenumber interpretation of the spectrum may be ap- 
propriate. An alternative view is that the fluctuations are not 
KAW-like f49l l69ll . in which case Taylor's hypothesis may 
break down. 

Background plasma parameters for the interval were deter- 
mined from the FGM Ir70ll and ESA l63ll instruments: so- 
lar wind speed i> sw = 320 km/s, magnetic field strength B 
= 5.5 nT, ion number density rii = 16 cm -3 , ion perpendic- 
ular temperature Tjj = 9.0 eV, electron perpendicular tem- 
perature Ti e = 11 eV, ion temperature anisotropy (T±/T\\)i 
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FIG. 3. Histogram of electron density spectral index values between 
ion and electron scales (3 < kpi < 15). 



= 0.90 and electron temperature anisotropy (T±/T\\) e = 1.0. 
The Doppler shifted kinetic scales calculated from these pa- 
rameters are marked in Fig. |2](under Taylor's hypothesis) with 
vertical dashed lines. 

In Fig. [2] at large scales (2xl0 -3 Hz - lxlO -1 Hz) a 
power law spectrum can be seen that is consistent with pre- 
vious measurements of the spectral index being around -5/3 
JH HHH. Just before the ion scales (0.1 Hz - 0.7 Hz) 
the spectrum flattens, which has also been seen previously 
1 55 , 5(1 71 ] and has been attributed to the turbulence becom- 
ing compressive Il4lll72ll or to pressure anisotropy instabilities 
17311 . The presence of these features suggests that the mea- 
surement technique is working well. 

After steepening at ion scales, the spectrum flattens to a 
constant value for / sc > 100 Hz. This is roughly consis- 
tent with the expected instrumental noise level and has been 
marked as a dotted line in Fig. |2] The spike near 1 kHz is of 
unknown origin but is not important for this analysis. Between 
ion and electron scales a power law can be seen. The spectral 
index over the range 3 < kp\ < 15, where the signal-to-noise 
ratio is large, was calculated from the gradient of the best fit 
line in log-log space and found to be -2.7. For greater ac- 
curacy, this spectral index was calculated from the spectrum 
in Fig. [2J5 and not the smoothed spectrum in Fig. [2};. Other 
spectral estimator techniques were used (e.g., windowed and 
wavelet transforms) with similar results. Since the spectrum 
reaches the noise floor around electron scales it is not yet pos- 
sible to determine whether it steepens, flattens or remains the 
same here. 

Spectral index variability. — The same procedure was ap- 
plied to 16 other intervals from October 2010 to January 201 1 
of duration between 6 min and 21 min. All of the intervals 
contained slow wind with 290 km/s < v sw < 350 km/s. Due 



to various sources of variability 11651 1741 . a different calibra- 
tion curve such as in Fig. Q~]was generated for each interval 
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FIG. 4. Spectral index (3 < kpi < 15) as a function of ion temper- 
ature anisotropy and ion parallel beta. The upper dashed line is the 
mirror instability threshold and the lower dashed line is the oblique 
firehose instability threshold from lT75ll . 



from a few hours of data containing the interval. A histogram 
of all 17 spectral indices is shown in Fig. [3] The mean spectral 
index is -2.75 ± 0.06, where the error is the standard error of 
the mean. 



It has been noted 15211 that studies of turbulence at kinetic 
scales in the solar wind should consider the contribution of 
instability generated fluctuations to the power spectrum. For 
example, the power in magnetic field fluctuations at the ion 
gyroscale is enhanced during times when the solar wind is 
marginally unstable to the firehose and mirror instabilities 
1 5211 . To examine their possible effect on the measurements in 
this Letter, the intervals were plotted in the instability param- 
eter space and colored according to the spectral index (Fig. 0). 
It can be seen that there is no consistent trend of spectral index 
with proximity to the thresholds, suggesting that the spectral 
indices measured here are not affected by these instabilities. 
A larger survey, with more coverage of the parameter space, 
however, would be required to make a more general statement. 

The flattening of the density spectrum above ion scales is 
present in all of the intervals reported here, irrespective of 
their location in Fig. [4] This suggests that it is inherent to 
the turbulent cascade, rather than being due to the pressure 
anisotropy instabilities. It it consistent with interpretations 
that the flattening is due to the compressive KAW fluctuations 
starting to dominate the density spectrum as the ion scales are 
reached H f72 1 . 

Discussion. — The spectral indices of density fluctuations 
measured here are similar to those obtained in measurements 
of the magnetic field at these scales (e.g., Jl, 15]). In partic- 
ular, the mean density spectral index of -2.75 ± 0.06 is the 
same to within errors as the universal magnetic field spectral 
index of -2.8 proposed in [15]. This is consistent with a cas- 
cade of fluctuations in which magnetic field and density are 
coupled and have the same spectral index, such as KAW tur- 
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The measured density spectrum, however, is steeper than 
the prediction of -7/3 for a pure whistler or KAW cascade. 
This has also been seen in 3D simulations of the magnetic field 
spectrum of both whistler [40Q and KAW IB6I1 turbulence that 
include kinetic effects and also a recent fluid simulation 13811 . 
There have been several explanations for the steep spectra, 
which rely on either energy being damped from the cascade 
or the intermittent nature of the fluctuations [ 26 34 - 3fl . In 
particular, the measurement is close to the specific spectral 
index of -8/3 predicted in l38ll . Other possible explanations 
include the applicability of Taylor's hypothesis (see earlier) 
and anisotropy of the scaling with respect to the mean field 
direction iflil 

Recently, the density spectrum in the Earth's foreshock re- 
gion has been measured at higher frequencies from 7.7 Hz to 
152 Hz and the perpendicular spectrum was reported to have 
a spectral index of -1.6 117611 . These results cannot be directly 
compared to the results of this Letter or to the dispersive cas- 
cade predictions, which are for solar wind turbulence above 
electron scales, and these shallow large amplitude spectra re- 
main to be explained but may be related to foreshock pro- 
cesses or the measurement technique. 

Finally, we note that compressible turbulence in general is 
poorly understood, even in neutral fluids. It has recently been 
proposed that compressibility would cause the energy trans- 
fer rate to vary locally H77I1 . which was suggested to explain 
results from compressible hydrodynamic turbulence simula- 
tions 17811 . We, therefore, have some way to go to fully un- 
derstand the nature of compressible plasma turbulence. The 
measurements in this Letter place an important constraint on 
theoretical descriptions of such turbulence and the calibration 
technique will allow the possibility for more detailed analyses 
of solar wind density fluctuations at these scales. 
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